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Abstract 

We introduce a new two-dimensional model with diagonal four spin exchange and an exactly 
known ground-state. Using variational ansatze and exact diagonalisation we calculate upper and 
lower bounds for the critical coupling of the model. Both for this model and for the Shastry- 
Sutherland model we study periodic systems up to system size 6x6. 
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I. INTRODUCTION 



In this paper we study two frustrated two-dimensional quantum spin models with exactly 
known singlet dimer ground states. 

First we consider the Shastry-Sutherland model Q , a two dimensional Heisenberg model 
with an additional frustrating coupling of strength Ji on every second diagonal bond (see 
Fig^) with Hamiltonian 



H = J2Y. ^^^1 + ^1 Yl ^^^r (1) 

<i,j> <i,j>diagonal 

Though the model shows similar properties also for higher spins, here we restrict ourselves to 
the case of spin i, and the Si denote spin-1/2 spin operators located at site i. The model Eq^ 
has attracted attention recently, when it was suggested, that the magnetic properties of the 
substance SrCu2 (603)2 [3] are very well described by its frustrated spin-spin interactions. 

The second model we investigate incorporates a diagonal ring exchange term on every 
second plaquette (see Fig|21). It is given by the Hamiltonian 



H = J2'^SiSj + Ji ^ {- + Sk^Sk3){- + Sk^Sk^,) (2) 

<i,j> plaquette k 

where again the Si denote spin-1/2 spin operators and the second sum runs over every second 
plaquette k as depicted in Fig|21 The four-spin interaction is normalized in such a manner 
that the dimer ground states, which play a major role in what follows, have zero energy. 
The extra four spin interaction of Eq|21 can be viewed as a "truncated" version of the ring 
exchange which has been discussed in the context of the two-dimensional mother substances 
of High a a„d_wKich have been found to pla. an essential .o.e in spin .adde. substances 



like Sri4Cu2404i 



Both models display rich phase diagrams and both have exactly known singlet dimer 
ground states. The latter property is otherwise rarely encountered when dealing with two- 
dimensional quantum spin systems: The ground state of the Shastry-Sutherland model is 
non degenerate and consists of (singlet) dimers located along the diagonal bonds, for the 
plaquette model the corresponding ground state is two-fold degenerate, and each ground 
state corresponds to a covering of the two types of diagonals with dimers. 



To facilitate the following discussion we define for both models the variable x := J2I J\i 
which can be viewed as an inverse frustration. The models undergo a zero temperature 
first order phase transition to the exactly known singlet dimer state at a critical inverse 
frustration x^- The numerical value of Xc is different for the two models. 

The aim of this paper is to introduce the basic properties of the plaquette model and 
to determine the value of Xc with best possible accuracy for both the Shastry-Sutherland 
and the plaquette model. Hereby we want to demonstrate, that also for two-dimensional 
frustrated models which are hard to tackle numerically, by use of the Lanczos technique 
on finite clusters not only reliable but also precise results for the infinite system can be 
obtained. We also hope, that our results are of use to gauge methods using uncontrolled 
approximations, whose errors are usually hard to assess. 

The paper is organized as follows. SectlUlis devoted to a study of the plaquette model 
using variational methods and finite cluster analysis. In particular two types of variational 
ansatze are used to give upper and lower bounds for the models critical coupling Xc- In 
Sect JIII A1 we present a study of the Shastry-Sutherland model of unprecedented precision 
using the Lanczos algorithm for system sizes up to 6 x 6 with periodic boundary conditions. 
In Sect JIII B1 we present a similar analysis for the plaquette model on periodic systems and 
give a comparison of the two models. In Sect JI VI we summarize our results and give an 
outlook to future possibilities. 

II. THE PLAQUETTE MODEL 

A. The ground state phase diagram 

Some first insight in the model can be gained by considering the smallest possible cluster, 
which consists of a plaquette of four spins only. This system has two singlet eigenstates 
with spin S" = 0, three triplets with = 1, and one "ferromagnetic" state with S = 2. 
The ground state diagram can be easily drawn for this system (see FigOl). Moving from 
the singlet dimer phase in anticlockwise direction at x = ^ = | one enters a second singlet 
phase usually referred to as antiferromagnetic phase or Neel phase and at J2 = 0, Ji < 
starts the ferromagnetic phase, which at x = — 1 meets again the singlet dimer phase. For 
iV = 4 the singlet dimer phase is degenerate with two S = \ phases. The third S = \ state 
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FIG. 1: Shastry-Sutherland model with J2 along diagonal bonds and Ji along vertical and hor- 
izontal bonds. The dimer ground state is build up by dimers i.e. "^(Ti ~ it) pairs along the 
diagonal bonds. 
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FIG. 2: Plaquette model with four spin interaction. 

is a ground state only for J2 = and Ji < 0, where it is degenerate with the Neel and the 
ferromagnetic state. 

It is still an open question to what extent a third phase phase between the dimer and 
the Neel phase exists for the infinite system. This problem is similar to the still unsolved 
issue in the Shastry-Sutherland model, for which various scenarios of intermittent phases 
between the Neel phase and the singlet dimer phase have been discussed. 

The phase boundaries between the Neel phase and the ferromagnetic phase and between 
the ferromagnetic and the dimer phase are exact, i.e. they coincide with the boundaries 
of the infinite system. However the value xi^~^^ = 0.5 of the phase boundary between the 
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FIG. 3: Ground state phase diagram of the plaquette model. The phase boundaries between the 
dimer and the Neel phase is drawn at J2 = Ji/2 which corresponds to a system with = 4 spins. 

dimer and its adjacent phase is not correct for the infinite lattice. To determine the value 
of Xc also for the infinite lattice, we calculated both lower and upper bounds and also study 
finite clusters with periodic boundary conditions. 

B. Lower and upper bounds for the Plaquette model 

Using variational ansatze in Ref. lower and upper bounds for the critical values Xc 
were discussed for the Shastry-Sutherland model. Here we perform a similar analysis for 
the plaquette model. As was pointed out by P.W. Anderson p] the ground state energy of 
open clusters Hi with no overlapping bonds, which cover the whole plane, represent a lower 
bound for the ground state energy of the infinite system. Similar to the Shastry-Sutherland 
system we can infer thus from the 4 spin plaquette that Xc^"* = 0.5 is a lower bound for 
Xc of the infinite system. An analysis of systems of larger size (see FigE]) gives as best 
lower bound Xc > 0.54263, which is the critical coupling of an open system of size = 31. 
Unfortunately the unsystematic variation of the bound with the system size does not allow 
a meaningful extrapolation. 

To obtain an upper bound on the critical coupling Xc the Hamiltonian is split into clusters 
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FIG. 4: Xc calculated on finite clusters with open boundary conditions for the Shastry-Sutherland 
model Eq^and the plaquette model Eq|21 The values for the Shasty-Sutherland model are taken 
from Ref.jsl and are listed here for comparison. 

without common spins and into external bonds connecting the clusters 

H = -f^dLster + Hhond- (3) 

clusteT 

This is graphically depicted in FigElfor clusters with four spins. As variational state we use 
a product state Hdnster l^c«), where \^ci) is the ground state of the cluster with Hamiltonian 
H duster- coutrast to the Shastry-Sutherland model however, where the expectation value 
of the external bonds vanishes for the ground state, which has spin 5 = 0, in case of the 
plaquette model the expectation value 

( n '^ciw.on^ n (4) 

duster duster 

has to be explicitly evaluated. Again a first bound can be found by considering clusters of 
size iV = 4. Best results are obtained here by using a plaquette with no four spin interaction 
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FIG. 5: A variational configuration of the "good type", consisting of two empty plaquettes with 
... . W.. „, . .p. ... i. P„ . 

corresponding configuration of the Shastry-Sutherland model |5|| one obtains an upper bound of 
|, wheras a configuration with the diagonal coupling Ji on the plaquettes ("bad type") gives only 
I as upper bound. 
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FIG. 6: Xc calculated on finite clusters with periodic boundary conditions for the Shastry- 
Sutherland model Eq^and the plaquette model EqEl 
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FIG. 7: Xc plotted versus for the Shastry-Sutherland model Eq^ (triangles and circles) 

and for the plaquette model Eq|21 (squares and circles). The lines connect systems with the same 
shape 

as basic cluster (as shown in FigEI). One thus obtains a^c < H ~ 0.6458. Again the 
numerical value can be improved here by considering larger clusters. Since the situation is 
more complicated in the case of the plaquette model, because one needs to consider also the 
bonds, we restrict ourselves here to a square lattice of size = 16, from which a numerical 
diagonalisation gives Xc < 0.627(1). 

To conclude, we find 0.54263 < Xc < 0.627(1) as exact bounds on the critical coupling of 
the plaquette model. 

III. THE SHASTRY SUTHERLAND AND THE PLAQUETTE MODEL ON SYS- 
TEMS WITH PERIODIC BOUNDARY CONDITIONS 

A. The Shastry Sutherland model 



We do not intend to give a detailed description of the Shastry-Sutherland model here since 
over the years it has been vastly discussed (see e.g. . In this paper we summarise only 

)riefly the results on Xc, which are most relevant for our work based on Lanczos technique 



Using variational approaches |5| it can be inferred from clusters of 31 and 32 sites, that 
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0.5914 < Xc < 0.6955. Thus by calculating variational bounds the margin left open for Xc is 
approximately 0.1. 

In this paper we discuss finite systems with periodic boundary conditions, which usually 
are closer to the infinite size system than systems with open boundaries. Also they allow to 
study larger systems, because translational invariance can be used to further reduce the size 
of the subspaces (for a list of the symmetries used for the different systems see table ITTT]) . 

Here we concentrate in particular on a precise determination of the critical coupling of a 
6x6 system. To successfully study this problem we had to deal with subspaces of dimension 
5 ■ 10^. We accomplished this by using parallel OpenMP programming techniques. We also 
investigated a = 32 square system and periodically closed 4 x n systems with up to n = 8. 
An overview of all the clusters with periodic boundaries is given in FiglHl Our results for 
the 6x6 system are shown in table |n] and for the A^ = 32 square systems in table HI 

As a technical point we mention that except for the trivial case of A^ = 4, there is no 
argument why the ground state energy for x > Xc should be linear in x. To obtain optimal 
results we have therefore chosen values of x as close as possible to Xc (see table n and HH). 
This is most important for the A^ = 36 system, for which calculating extra points is time 
consuming. We thus find as result for A^ = 36 a critical value of Xc = 0.674408, for the 32 
square we obtain Xc = 0.67724 and for the 4x8 stripe Xc = 0.66509. More values for smaller 
systems can be found in FigEl 

When studying two-dimensional clusters with periodic boundary conditions it is obvious, 
that the finite size effects depend not only on the number of spins in the system, but also 
on the geometry of the clusters. Here we have two types of systems, one consists of the 
A^ = 4 X n stripes with A^ = 16, 24, 32 spins and square systems with A^ = 16, 36. As 
was already pointed out in the stripe type systems with open boundaries show a linear 
finite size behaviour in 1/A^. We therefore extrapolate the A^ = 16,24,32 stripe systems as 
a function of 1/A^ and find Xc = 0.6624, which is an extrapolation for the 4 x n— stripe, 
however it cannot be taken as a limit for the two-dimensional system, because we deal 
essentially with a one- dimensional system. 

On first sight it seems a good idea to try an extrapolation of the A^ = 8, 18, 32 systems. 
This is however not feasible, because the A^ = 18 system, has rather pathological periodic 
boundary conditions and cannot be taken as a square system. 

The line connecting the system with A^ = 16 and A^ = 36 is not an extrapolation, since 
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the dependence on which is motivated by j8|,|9|,|l0| is a guess. The hne crosses the 

ordinate at Xc = 0.678, assuming a 1/N dependence gives Xc = 0.681 and assuming 
yields Xc = 0.6764. All three values lie within the upper and lower bonds from Ref.jsf. 

The fact, that the Xc for = 16 is smaller than for = 36 prompts us to believe that 
the true Xc should be larger than Xc = 0.6744. If true, this means, that the value Xc = 0.63 
fouud by Knetter et, a. Q the flow equafou method ,s too s,„a.L We roughly agree 
with the value Xc = 0.7±0.01 from Ref. |12|. which is obtained by extrapolating results from 
systems with 4 — 20 spins. 



B. The Plaquette Model 

The finite size behaviour of the plaquette model is in many respects similar to the Shastry- 
Sutherland model. Again the 4 x n stripes have lower Xc than the two-dimensional systems. 
Extrapolating the A^ = 18 and A^ = 32 squares in the above fashion gives Xc = 0.627 which 
coincides with the upper bound obtained from the A^ = 16 system. A 1/N dependence 
is not compatible with the upper bound and a 1/A^^ extrapolation gives Xc = 0.616. It 
seems more reliable to draw conclusions from the A^ = 36 and the A^ = 16 system. For the 
A^ = 6 X 6 square system we find Xc = 0.59136, which is again larger than the corresponding 
Xc = 0.56198 of the 4x4 system. We expect (though again there is no rigorous argument) 
that the critical coupling of the infinite size system should be larger than Xc = 0.59136. The 
line connecting the point of the 4x4 and 6x6 when plotted versus 

l/Ar3/2 

crosses the 

ordinate at Xc = 0.603 (assuming a 1/A^ dependence gives Xc = 0.615 and assuming 1/A^^ 
yields Xc = 0.599), all values lie between the rigorous upper and lower bonds of SectHTl 



IV. SUMMARY AND OUTLOOK 



We have introduced a plaquette model with an exactly known ground state and calculated 
the critical coupling Xc of the dimer phase in the model using the Lanczos algorithm for 
systems up to A^ = 36. For a 6 x 6 system we find Xc = 0.59136, as rigorous bounds for the 
critical coupling of the infinite system we obtain 0.54263 < Xc < 0.627(1). 

We also calculated the critical coupling of the Shastry-Sutherland model and obtained 
Xc = 0.6744 for a system with periodic boundaries and size A^ = 6 x 6. Extrapolating 
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J2/J1 
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J2/J1 
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0.67700000 


-0.374999999998987 
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-0.374999999999945 
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-0.374999999997117 


0.60000000 
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0.61000000 
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0.67724600 


-0.374999999997179 
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-0.375000015810957 






0.67730000 
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0.67750000 


-0.375115180019318 






0.67780000 


-0.375252606784488 






0.68000000 


-0.376271688712272 






0.70000000 


-0.386720856118170 







TABLE I: Data points calculated for the square = 32 Shastry-Sutherland model (left) and 
plaquette model (right). The dimension of the subspace containing the ground state is 75.164.451. 



J2/J1 


Eo/N 


J2/J1 


Eo/N 


0.50000000 


-0.374999999999425 


0.50000000 


8.16485767884956e-13 


0.66697792 


-0.374999999995975 


0.59135000 


5.54148393838716C-12 


0.67441000 


-0.375003966117052 


0.59140000 


-1.701097037248866-05 


0.67500000 


-0.375286177720292 


0.59330700 


-0.000846068815337997 


0.68000000 


-0.377704740424690 


0.59331000 


-0.000847382217467496 


0.70000000 


-0.388578449404367 


0.62000000 


-0.0133823277140735 


0.80000000 


-0.451753105004501 


0.63000000 


-0.0183954992316725 


0.85000000 


-0.484334145942023 







TABLE II: Data points obtained for the A^ = 6 x 6 = 36 Shastry-Sutherland model (left) and 
plaquette model (right). The dimension of the subspace with the ground state is 504.174.594. 

two-dimensional systems is subtle. Our best guess for the infinite system is Xc = 0.678 for 
the Shastry-Sutherland model and Xc = 0.603 for the plaquette model. 

We hope that our work will prompt further investigations of the plaquette model, in 
particular a detailed comparison of the magnetic properties of the plaquette model with 
the "full plaquette interaction" expected in high mother substances seems an attractive 
problem. 
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system 


dimension of subspace 




momentum 


reflection 


spin inversion 


32 square (periodic) 


75.164.451 


X 


X 




X 


4x8 (periodic) 


37.582.307 


X 


X 




X 


6x6 (periodic) 


504.174.594 


X 


X 




X 


28 (open) 


20.058.300 


X 






X 


31 (open) 


150.283.611 


X 




X 





TABLE III: Table of the dimensions of the subspaces and symmetry operations, which were im- 
plemented for the 32 square, the 4x8 and the 6x6 system with periodic boundaries and for the 
28 and 31 systems with open boundaries. 
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